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Electrophoretic mobility without charge driven by quadrupolar polarization of the nanoparticle/water interface 

Dmitry V. Matyushov^'0 

Center for Biological Physics, Arizona State University, PO Box 871504, Tempe, 
AZ 85287-1504 

Quadrupolar polarization of the interfacial water layer, typically found for water at contact with hydrophobic 
solutes, couples with the uniform external field to produce the force acting on a suspended particle. This force 
exists even in the absence of a net particle charge and its direction is consistent with an effectively negative 
surface charge. The surface charge density scales as with the particle radius R, and the corresponding 
contribution to the ^-potential scales as R~^. The charge density and C-potcntial are consistent with the 
values typically reported for oil drops and bubbles in water at i? ~ 3 — 10 nm, but decay below those values 
for larger particles. The study indicates that electrophoretic mobility without charge is possible, but is limited 
to relatively small particles a few nanometers in size. 
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I. INTRODUCTION 

The standard approach to electrokinctic mobility of 
particles suspended in solutioi>i starts with the assump- 
tion that the force acting on a particle carrying total 
charge q is given by the product of q and the Maxwell 
electric field E 



F = qE. 



(1) 



Mobility resulting from this force is a complex function of 
the surrounding electrolyto^. and hydrodynamic bound- 
ary conditions.— Significant literature is devoted to the 
subject,—!^ the current contribution is focused solely on 
deriving the force acting on a particle a few nanometers 
in size. 

Despite its simplicity, Eq. ([1]) carries a number of ap- 
proximations. First, in most practical situations the 
overall charge is distributed over the particle's surface. 
Replacing a generally nonuniform surface charge density 
with the total charge q is justified only if the surface 
of the particle coincides with the equipotential surface, 
which is, for instance, the case for a metal particle. The 
assumption of equipotential surface is much less satis- 
factory for insulating surfaces, especially at low charge 
densities. Therefore, even for a particle in vacuum, the 
force acting on it will deviate from Eq. ^ if the parti- 
cle is non-conductive and its surface charge is caused by 
adsorbed ions. 

An alternative to assuming the overall charge localized 
at a point, or distributed over an equipotential surface, 
is to directly calculate the force by integrating the stress 
tensor over the particle surfacci^ This approach connects 
the external field with the surface charge density, such 
that the cross term between the former and the latter 
leads to the dragging force (see below). While this formu- 
lation gives an expression identical to Eq. ([ij for a metal 
particle, the result generally differs from Eq. ([T|) for an in- 
sulating particle. The non- uniform distribution of surface 
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charge of an insulating particle can be expanded in terms 
of its orientational components (Legandre polynomials 
for axial symmetry). The first-order expansion leads to 
the dipolar polarizability of the particle and its interface, 
describing the dipole induced at the particle in response 
to an external field and including both electronic and 
permanent-charge susceptibilities. This induced dipole 
is not bound to the assumption of equipotential surface, 
implicit to Eq. ([1]), and thus allows a tangential com- 
ponent of the field at the particle surface. Importantly, 
the polarizability attached to the solute accommodates 
not only the dipole moment intrinsic to the solute, but 
also the interfacial dipolar polarization, which always ex- 
ists at a dielectric interface^ and has recently received 
attention in connection to the problem of electrokinctic 
mobilityi^r— 

Once the idea of surface charge density is adopted for 
the calculation of the dragging force, it opens the poten- 
tial for introducing molecular properties of the interface 
affecting the overall surface charge density. Specifically, 
polar liquids in general and water in particular show sig- 
nificant changes of their local density (wetting or dewet- 
ting) at interfacesii and, in addition, specific patterns 
of orientational order J^"— It was found that interfacial 
dipoles tend to orient parallel to the interfacial plane 
both for simple dipolar fluidsi^"— and for wate r^'^'^^i^^ 
at contact with a nonpolar wall. This peculiar orienta- 
tional structure is also found at the air-water interface 
by computer simulations^^ and, experimentally, by sec- 
ond harmonic generation . ^"^i^^" — 

Here we calculate the force acting on a nanometer par- 
ticle in a uniform external field by surface integration of 
the electrostatic (Maxwell^i^) stress tensor. The main 
question addressed here is whether this more precise ap- 
proach can offer significant changes to the standard ex- 
pression for a point charge given by Eq. ([1]). 

We start with introducing the dipolar polarizability of 
the particle into the equations for the force. This addi- 
tion leads to the equation for the force distinct from Eq. 
(H]). While the standard result is obtained for a metal 
particle, the dragging force generally depends on the in- 
terfacial solvent structure through the polarizing cavity 
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field. Depending on the structure of the interface, either 
the Maxwell scenario, screening the solute polarizability, 
or the Lorentz scenario, elevating the effect of the solute 
polarizability, takes placci^ Altering the solute polariz- 
ability, by for instance photoexciting electron-hole pairs, 
can potentially discriminate between the possible scenar- 
ios. 

We next investigate the role of the interfacial 
quadrupolar polarization and show that it projects itself 
on a second-order, in Legendre polynomials, component 
of the surface charge density. This charge density, unre- 
lated to the density of trapped surface charges, couples 
in the Maxwell stress tensor to the external field, thus 
yielding an additional component of the force acting on 
the particle. The main result is that the external field 
couples to the effective charge 

= (7e + §5(72 (2) 

to be substituted into Eq. ([T]) in place of the Coulomb 
charge q (with an additional correction originating from 
the particle polarizability, as discussed below). The effec- 
tive particle charge in Eq. Q is modified compared to q is 
two essential ways. First, q is replaced with the electroki- 
netic charge q^ which includes the compensating charge 
of the electrolyte within the shear surface enveloping the 
part of the liquid moving together with the particle.^^ 
Second, the second summand introduces an additional 
surface charge density a2 arising from quadrupolar po- 
larization of the interfacial water; S is the surface area. 
Quadrupolar polarization appears spontaneously in the 
interfacial layer driven by the necessity to accommodate 
both the dipole and quadrupole moments of the inter- 
facial waters to minimize the free energy density of the 
interface.—"— As we show below, this spontaneous polar- 
ization, existing even in the absence of external polarizing 
fields, can produce an electrostatic potential and associ- 
ated with it surface charge density. The effective charge 
is non-zero and, as we show below, negative at g = 0. 
This result is the presence of a dragging force acting on 
a particle of zero charge and caused by the interfacial 
structure. 



II. NANOPARTICLE IN A DIELECTRIC 

In order to set the stage for the theory development, 
we will start with the simple case of a spherical particle 
immersed in a dielectric with the dielectric constant e. 
The nanoparticle is assumed to carry the uniform surface 
charge density ctq and to possess the dipolar polarizability 
OfQ. The entire system, composed of the nanoparticle 
and the surrounding dielectric, is placed in an external 
uniform electric field Eq aligned with the z-axis of the 
laboratory coordinate frame (Fig. [T]). 

To illustrate our derivation steps, we first simplify the 
problem even further by removing the dielectric and plac- 
ing the particle in vacuum. We therefore start with e — 1. 
The total force acting on the particle along the z-axis is 




FIG. 1. Cartoon of the water molecule at the surface of a 
spherical particle. The normal to the surface forms the angle 
6 with the direction of the external field along the z-axis of 
the laboratory coordinate frame. The orientation of the wa- 
ter molecule is specified by two angles: 9^ (between the sur- 
face normal and water's dipole moment) and x-w (between the 
plane containing n and n and the plane of the water molecule). 

obtained by surface integration of the Maxwell stress ten- 
sor (Jik, contracted with the Cartesian components fii of 
the outward surface normal vector— 

F, = a,,n,dS, (3) 

where dS is the surface area differential and summation 
runs over the common Cartesian indexes. 

The electric field E entering the Maxwell stress tensor 
CTj^ = (47r)"^ (EiEk - {l/2)SikE'^) is calculated from the 
electrostatic potential (f>{r) outside the spherical particle 

<^(r) = -^o^cos0(l-^) (4) 

in which q = AirR^ao is the overall charge and R is the 
radius. Further, is the angle between the normal fi, 
outward to the surface of the particle, and the z-axis 
associated with the external electric field (Fig. [T]). 

The standard solution of the Maxwell boundary prob- 
lem for a metal particle^ results in ao — R"^ in Eq. ([4]). 
The electric field is then normal to the surface at r = i?. 
Movements of charges trapped at the surface of an insu- 
lating nanoparticle do not necessarily follow the rules of 
a metal conductor, and so the dipolar polarizability ao 
accounts for possible deviations from the charge distri- 
bution characteristic of a metal particle. In particular, 
the electric field has a tangential component at the par- 
ticle's surface when ao ^ R^ . Note also that ao does not 
necessarily represent the electronic dipolar polarizability 
of the particle, but can also account for the distribution 
of permanent charge. It therefore should be associated 
with the Debye polarizability,— or dipolar susceptibil- 
ity, connecting the total dipole induced at the particle 
to the external electric field Eq- This overall dipole in- 
cludes both the polarizable electrons and the permanent 
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charges within the particle responsible for a permanent 
dipole aligned along the field. 

In order to clarify the origin of the force acting on a 
finite-size particle, we apply the connection between the 
normal component of the field at the surface En and the 
surface charge density, Aira = En- This relation leads to 
the following expression for the surface charge density 



a{9) = (To + T-£^o cos( 
47r 



2aQ 

i?3 



1 + ^ (5) 



In addition, one can re- write the expression for the force 
in terms of the normal component of the field only, with 
the resulli^ 



=Att^R^ j a{ef cos 9 d COB e. 



(6) 



The only non-zero term in this equation is the cross term 
between the uniform charge density of the trapped sur- 
face charges ctq and the polarization induced by the ex- 
ternal field, (X cosO. Integration in Eq. ([6]) yields 



F, 



(7) 



The standard result for the force acting on a point 
charge, F^ = qEq, follows when the surface of the particle 
coincides with the equipotential surface, i.e., when the 
particle polarizability corresponds to the limit of a metal 
sphere, ao = R^. If the particle is non-polarizable, uq = 
0, one gets only one-third of the force, Fz — {q/3)Eo. The 
range of possible scenarios broadens further when a polar 
solvent with the dielectric constant e > 1 is introduced. 

The simplified case of a particle in vacuum is presented 
here to stress that the standard Coulomb force acting on 
a uniformly charged sphere arises from the cross term 
between the surface charge density and the dipolar (first 
Legcndre polynomial) polarization of the particle. We 
will use this result further below to introduce the force 
in the absence of the overall charge, which arises from the 
cross term between the dipolar polarization of the parti- 
cle and an additional surface charge density arising from 
quadrupolar polarization of the hydration layer. The lat- 
ter is present even in the absence of the external field be- 
cause of the preferential orientational structure of waters 
in contact with the solute. 

Assume next that the spherical particle with the sur- 
face charge density (Tq and polarizability ao is placed in 
the dielectric with the dielectric constant e. The electro- 
static potential in Eq. ([4]) changes to 



-Er cos 6 1 



Xc 3 



_9 

er 



(8) 



where E = Eq/e is the Maxwell electric field in the bulk. 
Further, since ao represents the dipole induced at the 
particle in vacuum, a correction is required when the 



particle is placed in a dielectric. This is achieved by the 
cavity field susceptibility2ii2£ = Ec/Eq connecting the 
field inside the solute Ec ( "cavity field" ) with the field of 
the external charges Eq. 

When the dielectric constant eo can be assigned to the 
material of the solute, the standard boundary conditions 
of Maxwell's electrostatics require 



"0 



2e + eo 



(9) 



The dipole moment at the solute is therefore directed 
opposite to the external field when the solute is less po- 
larizable than the solvent (cq < e) and is along the ex- 
ternal field when the solute is more polarizable (eo > e). 
One again gets ao = R^ for a metal particle in the limit 
eo oo. Alternatively, if an empty cavity is introduced 
into the dielectric (eo = 1), one gets the standard ex- 
pression for the interface dipole moment Mq"' = aoE 
produced by polarizing the cavity interface^ 



3»o J 
"2e+ 1 



(10) 



where Qq is the cavity volume and P is the dipolar po- 
larization of the bulk. 

Equations ^ and ^ arc useful for connecting the 
present formalism to studies of macroscopic suspensions. 
However, the dielectric constant is not straightforward to 
define for nanoscale objects and the language of solute 
polarizability, commonly adopted for molecular solutes, 
is more preferable in such cases. 

Repeating the calculation of the overall force as shown 
above, one gets for the force acting on the particle 



1 



Xc 



2ao 
i?3 



(11) 



The contribution of the solute polarizability to the force 
is strongly affected by the cavity-field susceptibility. If 
the standard Maxwell prescription is used for this func- 
tion, Xc = 3eo/(2e -|- eo), the polarizability term in Eq. 

is strongly reduced by the screening introduced by 
a polar liquid, such as water. 

We have recently found that a different limit for the 
cavity field, designated as the Lorentz field, is reached 
for interfaces of polar liquids with nonpolar solutes^^i^ 
Xc = (e + 2eo)/(3e). The distinction between this out- 
come and the Maxwell limit is rooted in the difference 
between the orientational structure actually realized at 
interfaces of dense polar liquids and the one assumed in 
the dielectric models. While nonspecific freedom of the 
surface dipoles to reorient is assumed by the dielectric 
models, polar molecules, and water in particular, assume 
preferential in-plane orientations of their dipoles when 
in contact with a nonpolar surface*^ This orientational 
structure strongly disfavors dipolar response normal to 
the surface, thus suppressing the creation of the normal 
projection P„ of the dipolar polarization P at the di- 
viding surface 1^ The suppression of the normal polar 
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response alters the overall electric field inside the solute 
compared to the standard prescriptions. The solute po- 
larizability has a much stronger impact on the overall 
force acting on a charged particle in this latter scenario. 

The interfacial polar response, and thus the normal 
component of the interfacial dipolar polarization P„, 
is affected by the solvent structure formed around po- 
lar/ionic groups at the solute surfaccj^ Therefore, the 
actual polarization of the interface in real solutions will 
deviate from either of the two scenarios outlined above. 
Which scenario is realized can be established by altering 
the solute polarizability. This opportunity is particu- 
larly attractive for semiconductor nanoparticles. Their 
photoexcitation creates highly polarizable electron-hole 
pairs. The polarizability of an electron-hole pair can 
reach ^ 10^ A^^, and it scales as cx with the solute 
radius According to Eq. pT|) . the mobility of pho- 
toexcitcd nanoparticles should increase compared to the 
ground state, and the extent of enhancement gives direct 
access to Xc- 



III. QUADRUPOLAR POLARIZATION OF 
THE INTERFACE 

The above discussion highlights the general phe- 
nomenon of preferential orientation of surface molecules 
when two phases are in contact. Since water carries both 
a large dipole and a large quadrupole moments, a specific 
orientational order at the water interface gives rise to 
distinct dipolar and quadrupolar responsesi^ Those will 
have observable electrostatic signatures r^i^ and will also 
produce mechanical effects when movement of charges 
by the electric field is concerned. Here we show that 
the quadrupolar component of the interfacial polariza- 
tion couples to the external field to produce dragging 
force applied to a nanoparticle, which does not require a 
net electric charge. 

As mentioned above, only the normal component of 
the interfacial polarization creates surface charge density 
a (9). When the dipolar polarization of the interface is 
concerned, one therefore needs to consider the first-order 
orientational parameter in the surface layer pi = {Pi{ii ■ 
in)), where Pf(n • ih) is the Legendre polynomial as a 
function of the cosine formed by the normal to the solute 
surface n at the location of the water molecule and its 
dipole moment m; hats denote unit vectors (Fig.[T]). The 
statistical average here is typically performed over the 
waters residing in a few hydration layers closest to the 
solute. 

The order parameter pi defines the overall non- 
compensated radial dipolar polarization at the interface, 
and it is typically smalL ^^'^^i^^ While some residual po- 
larization might still exist, it will only contribute to the 
effective dipolar polarizability of the interface, which has 
been already included into the solute polarizability ag. 
This component will therefore be neglected for the re- 
mainder of our analysis. In contrast to pi, the second or- 



der parameter p2 = (f2(n ■ lii)) is non-zero and negative 
for water interfacing both nonpolar solutesSi and nonpo- 
lar planar surfaces This order parameter appears in the 
electrostatic potential created by surface quadrupoles, 
and one might anticipate that quadrupolar polarization 
is an important contributor to the overall polar response 
of the interface. 

Quadrupolar interfacial polarization adds a component 
proportional to the second-order Legendre polynomial to 
the surface charge density. The expansion of a {6) in Leg- 
endre polynomials then includes three terms 

2 

(T{e) =^aePe{cos0). (12) 

This equation can now be used in Eq. ([6]) to calculate the 
overall force projected on the external field 

F. = ^^(.„ + -.2j. (13) 
With cTi = {E/Ati){1 + 2xcao/R^), one gets 

where the effective charge of the solute is given by Eq. 

The main result of this derivation compared to Eq. 
(fTTj) is that the cross term between the first-order (dipo- 
lar) and second-order (quadrupolar) components of the 
surface charge density results in a non- vanishing drag on 
the particle. A non-zero force thus exists even at zero 
charge, q = 0. The model thus requires a distinction be- 
tween the isoelectric point, = 0, and the point of zero 
charge, q ~ 0. 

The derivation presented in the Appendix shows that 
the second-order expansion coefficient in the surface 
charge density in Eq. (|12p can be obtained from the 
quadrupolar polarization of the hydration shell as follows 

CT2 = imsh/R'') [p2Qzz + q2AQ] . (15) 

In this equation, Ngh is the number of waters in the hy- 
dration shell. The derivation in the Appendix also sug- 
gests that the radius R in this equation should be iden- 
tified with the radius of the surface including the hydra- 
tion layer producing a2- It is natural, in application to 
electrokinetic problems, to identify this layer with waters 
within the shear surface and then Ngh with the number of 
waters within the share surface. Since in most practical 
applications, the size of the suspended particle is signif- 
icantly larger than the diameter of the water molecule 
CTu,, we will make no distinction between the radius of 
the shear surface and the radius on the particle itself ap- 
pearing in the polarizability term. Given this association 
of the particle size with its shear surface, the electroki- 
netic charge should be assigned to the particle, as is done 
in Eq. ©. 
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Further, in Eq. ([15]), Qzz and AQ = Qxx — Qyy are 
the Cartesian components of the water quadrupole in 
the frame of molecular principle axes diagonalizing the 
quadrupole moment matrix. The water quadrupole is 
mostly non- axial, with AQ = 5.13 DxA and = 0.13 
D X A.— Since p2 and q2 in Eq. (fTC|) are of the same order 
of magnitude for nanomctcr-sizc solutes?^ the non-axial 
quadrupole component AQ dominates in Eq. (jl5p . 

Orientational order parameter p2 and q2 in Eq. (|15p 
are non-zero due to spontaneous orientational structure 
of the interfacial waters existing even in the absence of 
external fields. The order parameter p2 was mentioned 
above and is given by the average of the second Legendre 
polynomial of the angle 9^ between the surface normal 
and the water's dipole moment (Fig. [Ij. The second or- 
der parameter q2 oc (l?02 + ^0-2) in Eq. (|15p arises from 
the real part of the second-order rotational matrix -D^j„j 
transforming the quadrupole moment from the molecular 
to the laboratory frame of coordinatcsi^ It is therefore 
sensitive not only to the orientation of the water dipole 
relative to the surface normal, but also to the Euler angle 
Xw specifying the orientation of the water plane relative 
to the plane of the dipole moment and surface normal; 
XiD = when the two planes coincide (Fig. [T|). The re- 
sulting expression for this order parameter is (see the 
Appendix) 

(?2 = 5 (sin^ Oy, cos 2xw) ■ (16) 

The actual values of p2 and (?2 for hydrophobic and 
water-air surfaces point to in-plane alignment of the 
dipole moments.— The parameter q2 is additionally sen- 
sitive to the presence of dangling bonds, which occur 
about every fourth surface water facing a hydrophobic 
solute The typical conditions of hydrophobic solvation 
then suggest 92 < 0, while (72 — —0.15 is found in numer- 
ical simulations of Lennard- Jones solutes of about one 
nanometer in radiuSfSS This observation implies that (T2 
contributes to an effective negative charge of the solute 
when measured by electrophoretic mobility. 

IV. DISCUSSION 

The insertion of a particle in water creates specific ori- 
entation of waters at the interface. Because of molecular 
asymmetry, waters residing in the interfacial regions at 
opposite sides of the particle, along the direction z of the 
external field, respond asymmetrically to the field. The 
result is that the volume density of the electrostatic en- 
ergy, c>c -B(r)^, is different on the opposite sides of the 
particle, thus creating a gradient of the chemical poten- 
tial. This chemical potential gradient results in a non- 
compensated pressure, similar to the osmotic pressure 
arising from the gradient of chemical potential between 
two parts of the solution. The Maxwell tensor used in 
the present derivation replaces the free energy density 
with the surface stress^ which is easier to accommodate 
in specific calculations. 




water 



z 

FIG. 2. Cartoon of the plane-capacitor experiment in which 
the external electric field is created by two metal plates. The 
void in water, indicating the suspended particle, develops sur- 
face charges at the dividing surface separating the void from 
water. An excess of negative charge is meant to stress that an 
effectively negative charge density is present at the void-water 
interface. The interfacial charges will cause the correspond- 
ing image charges in the metal plates, thus altering the field 
produced by the capacitor. Shifting the void by the distance 
Az changes the overall distribution of the surface charge, thus 
altering the overall polarization free energy of the capacitor. 
This change in the free energy should be equal to the work 
done to move the particle. 

One still wonders about the thermodynamic balance 
when an uncharged particle is dragged by the force and 
finite work is done. The arguments here might be subtle 
and should involve the consideration of the electrostatic 
problem in a particular experiment, for instance in the 
plane capacitor experiment shown in Fig. [2l It is gener- 
ally accepted that dielectric interfaces induce interfacial 
charges. The creation of a void in the water dielectric 
creates surface charge density, which integrates into the 
overall surface dipole in Eq. (fTO|l .— It is also easy to real- 
ize that charges at the void's surface will produce image 
charges in the conductor plates, thus resulting in an alter- 
ation of the electric field inside the capacitor. The solu- 
tion of the corresponding Poisson problem will depend on 
the position of the void, and the entire free energy of the 
electrostatic polarization of the dielectric will be altered 
by the displacement of the void (Fig. [2|). This change 
of the electrostatic free energy goes into work done to 
drag the particle and is supplied by the external power 
source maintaining the constant voltage at the capacitor 
plates. Since ions do not screen dipole and higher inter- 
face multipoles, these arguments arc not affected by the 
Debye-Huckel electrolyte. 

Several parameters entering Eqs. (jl4[) - (|15p require bet- 
ter grasp of the electrostatic interfacial properties, which 
is mostly missing at the moment. To start off, the num- 
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ber of particle A^sh in the hydration sheh in Eq. p5l) 
is an empirical parameter when laboratory experiments 
are concerned, but can be calculated from atomistic nu- 
merical simulations. The size of the shell is determined 
by the depth of orientational order specific to interfa- 
cial waters. This penetration depth can be determined 
from the order parameter q2{r) calculated as a function 
of the thickness of the hydration layer r surrounding the 
nanoparticle. The saturation of q2ir) to its bulk value 
at r ~ (5 yields both the depth i5 of the orientational in- 
terfacial structure and, when the corresponding density 
pair distribution function is integrated over that region, 
the number of waters A^sh in the hydration shell. Assum- 
ing that the size of the nanoparticle is much larger than 
the depth of the hydration layer, one has A^sh — ^irR^pS, 
where p is an effective number density of the interfacial 
water. One can therefore rc-write the expression for ct2 
as follows 



0-2 ^ 7277^ 



5 92 AQ 
(Jw R"^ al, ' 



(17) 



where rjw = (7r/6)pcr^ ~ 0.41 is the packing density of 
water. With = 2.87 A and one gets an estimate for the 
effective surface charge density arising from quadrupolar 
polarization 



(72 ^ nie/R^)iq2S/a^). 



(18) 



For instance at q2 ~ —0.15^ and S/a.w = 3 one gets (72 = 
— e/nm^, where the radius of the particle is in nm. 
The value of —0.3 e/mn^ typically reported for interfaces 
between air bubbles and oil drops with wate r"^^'"^^ is thus 
achieved at i? ~ 3 nm and decays further with increasing 
the particle size. 

Experiment does not report the force, but, instead, the 
mobility of a particle given by the ratio of the station- 
ary velocity ve to the external fieldji p — ve/E. The 
mobility is linked to the (^-potential C by Smoluchowski's 
equation in the limit of large solutes, exceeding in size the 
Debye-Hiickel length of the electrolyte,— p, = Ceoff/(47r?7), 
where 77 is the viscosity (a factor of 2/3 appears in this 
equation if the size of the solute is below the Debye- 
Hiickel length^), ^-potential is the electrostatic poten- 
tial at the shear surface, i.e., the surface surrounding the 
solute at which the velocity of water starts to deviate 
from the velocity of the solute; all waters are assumed 
to move with the solute inside the shear surface (assum- 
ing wetting^) . Since these are typically waters with the 
energy of binding to the solute exceeding their thermal 
energy, an effective dielectric constant within the shear 
surface Ceft is used in the Smoluchowski equation)^ The 
concept of the local dielectric constant is not well defined 
in this context, not only because of conceptual difficulties 
of defining local polar response, but also because fluids 
become birefringent in the vicinity of a dielectric inter- 
face, i.e., polar response is anisotropic 

The mobility of colloidal particles is a complex function 
of the electrolyte concentration for a non-zero surface 



charge ctq 7^ 0. Our goal here is to estimate the magni- 
tude of C-potential produced by the interfacial quadrupo- 
lar polarization. It suffices to put ctq = for this purpose. 
The effect of the electrolyte can then be neglected and 
mobility determined by balancing the drag force with the 
viscous friction force GTrryugi?. We find for the zero-point 
charge (zpc) ^-potential 



_ IGnR 

Czpc — TT" — ^2, 



45e 



cff 



(19) 



where we have neglected the polarizability effect under 
the assumption of the Maxwell form for the cavity field; 
this gives the lower-boundary estimate for |Czpc|- Substi- 
tuting Eq. (HH) into Eq. (HH), one gets 



zpc 



-900 



(7„,ecffi?(nm) 



mV. 



(20) 



The parameters in this equation are not well established. 
However, (5 ~ 2 — 3 nm is typically accepted 1^— The 
value of —Czpc is then a few tens of mVs for a particle 
~ 5 — 10 nm in size. 



V. SUMMARY 

The standard equations for mobility in a uniform ex- 
ternal field require the particle to carry a net charge [Eq. 
(P)]. Dipole requires field gradient to produce the force. 
The same statement applies to the interfacial dipolar po- 
larization, either induced by the external field or existing 
in the interface in the absence of the field. We find no 
force arising from this source when the external field is 
uniform, as discussed previously in the literature.— How- 
ever, the force drag does appear from the coupling of 
a uniform external field to the quadrupolar polarization 
of the interface. This outcome is particularly important 
for aqueous interfaces since the symmetry of preferen- 
tial orientations of waters in nanometer-scale interfaces 
favors a substantial quadrupolar polarization ; ^'^'^'^ even 
in the absence of external polarizing fields. The present 
formulation thus predicts a non-zero electrokinetic effect 
for a particle carrying no charge. The direction of the 
force is consistent with an effective negative charge when 
the standard analysis of the electrophoretic mobility is 
applied. The C-potential produced by compensating the 
drag with hydrodynamic friction scales as R~^ with the 
particle radius R. It therefore become negligible when 
the particle size approaches (sub) micrometers and the 
traditional mechanism of charging the particle through 
accumulation of surface charge^i^ is expected to domi- 
nate in that range of sizes. 
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Appendix A: Derivation of Eq. lfT5|) 

We want to derive the electrostatic potential of charges 
Qi with coordinates Vi belonging to a molecule with the 
center-of-mass coordinates Tq located at the surface of a 
spherical solute. The geometry relevant to the problem 
at hand is shown in Fig. [3l and we assume that the po- 
tential is calculated at point r outside the solute, such 
that r > ro- The standard Coulomb expression for the 
electrostatic potential can then be multipole-expanded in 
powers of r^/r and ro/r. Since the expansion satisfies a 
number of rotational invariance constraints, the expan- 
sion in rotational invariants applies here^S 



E 



water 



r - ro - Ti 

i 



(Al) 

where Ygmi^) is the spherical harmonic, 
C{£i£2^;mim2m) is the Clebsch-Gordan coefficient, 
hats denote the unit vectors, and 



1 



21 + 1 



(47r)3(2£+ 1)! 



(24 + l)!(24 + l)! 



1/2 



(A2) 



In addition, £ = i?i -I- ^2 is imposed in the sum above to 
produce the right dimension for the electrostatic poten- 
tial. 

The sum over the molecular charges leads to the mul- 
tipolar moment expressed in the spherical coordinates 



Qim = ^ QiriYhnih)- 



(A3) 



Given that we are interested in the quadrupolar polar- 
ization of the interface, we put ^2 = 2 in Eq. (|Aip , with 
the result 



Ci£ 2 £ + 2;mim2m) 

Q2m2Ylmi (fo)^f*f2,m(r) 



(A4) 



The axial symmetry of the problem (Fig. [3]) suggests 
that the result should be invariant to rotations about z, 
thus requiring m = mi ~ 7712 = in Eq. (jA4p . The 
equation simplifies to 

<^2(r) = E 4§^?2o^«(fo)y;+2,o(f) (A5) 



with 



(£ + l)(^ + 2) 



(47r) 



5(2£+l)(2^ + 5) 



1/2 



(A6) 




FIG. 3. Cartoon of the solute (large sphere) with a water 
molecule (small sphere) at its surface, qi denotes an internal 
charge of the water molecule, is its coordinate relative to 
the water's center of mass, z-axis of the laboratory coordinate 
frame is also specified. 



We now take the average over the positions of the 
molecular quadrupoles at the surface of the sphere and 
their orientations in the absence of the external field. The 
spherical symmetry of the problem in the absence of the 
field requires that all positions at the surface of a spher- 
ical particle are equivalent, i.e., £ = 0. Next, in order to 
perform the average over the orientations of the molecules 
at the surface, we transform the quadrupole moment ma- 
trix from the laboratory frame to the molecular principle 
axes diagonalizing the molecular quadrupole. 

The system has spherical symmetry in the absence of 
the external field, and the orientation of the laboratory 
frame can be arbitrary. It is identified here by aligning 
its z-axis with Vq and its a;-axis with the plane formed by 
ro and the molecule's dipole moment m. The transfor- 
mation between the two frames is accomplished by the 
rotational matrix D^„{n) such that 



Q20-E^0n(^^)Q2n. 



(A7) 



Here, VL = (pwOwXw is a set of Euler angles that carry the 
rotation of the molecular frame into coincidence with the 
laboratory frame. Since to = in -D^„ in Eq. ([9]), <f)w = ^ 
and only two remaining angles perform the transforma- 
tion. The Euler angle 0^ gives the rotation of the water 
molecular dipole relative to the surface normal, while Xw 
is the subsequent rotation around the dipole moment m, 
which moves the water plane out from coincidence with 
the plane formed by the normal and m. 

Since the quadrupole moment in the rhs of Eq. (jA7p is 
defined relative to the principle axes of the molecule. 



Q20 —\ -;—Qzz, Q21 — 0, 
47r 



Q22 —Q2-2 




(A8) 



{Qxx Qyy) ■ 



Wc therefore obtain, with the account for the properties 
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of rotational matrices^^ (AQ = Qxx — Qyy) 



DLmQ2n 



where 



P2 = (P2(C0S6'^)), 

(72 = ^(sin^ 6'.„cos2x« 



gsAQ] , (A9) 



(AlO) 



and Pi {x) is the Legendrc polynomial. The angles 0^ and 
Xw are calculated relative to the normal to the surface at 
the position of the molecule (Fig. [T]). This is indeed how 
the orientational order parameters are calculated in the 
analysis of MD trajectories!^ 

Finally, one arrives at the electrostatic potential 



02 (r) 



„-3 



(All) 



where cos 6 = z ■ Vq and z is the unit vector specifying 
the z-axis of the laboratory frame (Fig. [3]). The surface 
charge density is obtained by multiplying this expression 
by the number of particles in the hydration layer and 
taking the radial derivative, 



(72(0)^2 (cos 0) = -N,i,dMr)/dr\ 



(A12) 



The derivation presented here also makes it clear that 
the potential 02 (r) is zero when r < tq. The surface at 
which (72 is calculated should therefore lie outside of the 
hydration layer producing the quadrupolar polarization; 
identifying R with the radius of the shear surface seems 
the right choice for electrokinetic problems. 
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